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In this paper, we investigate the dynamics of the Bloch vector of a single two-level atom which 
interacts with a single quantized electromagnetic field mode according to the Jaynes-Cummings 
model, where the field is initially prepared in a thermal state. The time evolution of the Bloch 
vector S(t) seems to be in complete disorder because of the thermal distribution of the initial state 
of the field. Both the norm and the direction of S(t) oscillate hard and their periods seem infinite. 
We observe that the trajectory of the time evolution of S(t) in the two- or three-dimensional space 
does not form a closed path. To remove the fast frequency oscillation from the trajectory, we take 
the time-average of the Bloch vector S(t). We examine the histogram of {S z (nAt)\n = 0, 1, N} 
for small At and large N. It represents an absolute value of a derivative of the inverse function of 
S z (t). (When the inverse function of y = S z (t) is a multi- valued function, the histogram represents 
a summation of the absolute values of its derivatives at points whose real parts are equal to y on the 
Riemann surface.) We examine the dependence of the variance of the histogram on the temperature 
of the field. We estimate the lower bound of the entanglement between the atom and the field. 
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I. INTRODUCTION 



The Jaynes-Cummings model (JCM) is a solvable quantum mechanical model of a single two-level atom in a single 
electromagnetic field mode [H, 0, S 0, S] • This model is originally designed for studying a spontaneous emission. The 
interaction term is obtained by the rotating wave approximation. In this interaction, each photon creation causes an 
atomic de-excitation and each photon annihilation causes an atomic excitation. 

If the photon number is sharply defined in the initial state of the field, the JCM shows the Rabi oscillations in the 
populations of the atomic levels. If the initial state of the field is a coherent state, the oscillation of the mean photon 
number collapses and revives in the JCM. In this way, the JCM reveals the quantum natures of the radiation. 

Because the JCM is exactly solvable, it is investigated from various viewpoints. The JCM whose boson field is 
prepared in a thermal state is discussed in Refs. [1,0, 9- Thermodynamics of the JCM is discussed in Ref. @. In 
■ this analysis, the grand partition function of both the atom and the boson field is considered. Extended JCMs are 
studied as dissipative models [H, [ll| . 

Recently, the JCM has been used for describing the evolution of the entanglement between the atom and the 
field [H], [l3| . In Refs. [H, EH , the electromagnetic field is assumed to be initially prepared in a thermal state. In 
these papers, the JCM is regarded as a source of the entanglement between the atom and the field. In Ref. [3], 
this idea is advanced and generation of the entanglement between two atoms interacting with a single-mode thermal 
field according to the JCM is discussed. In this model, the atoms which are initially in a separable state obtain the 
entanglement through the time evolution. 

In Ref. 15], evolution of the entanglement between an atom and a single-mode field described with the JCM under 
phase damping is studied. In Ref. [161 ] . an atom interacting with two cavity modes is considered and it is shown that 
this system can be reduced to the JCM. Using this system, generation of entangled coherent states is discussed. In 
Refs. [13, E], the entanglement in an extended JCM where two atoms interacts with a single-mode field is studied. 

As an interesting phenomenon related with the time evolution of the entanglement in the JCM, so called sudden 
death effect is studied [l9|, [H,[2l[. In these references, two isolated atoms, each of which is located in its own Jaynes- 
Cummings cavity, are assumed. If these atoms are in a certain entangled state initially, the entanglement disappears 
in a finite time. This phenomenon is experimentally demonstrated in Ref. [22| . In Ref. [23j . an attempt to find 
invariant entanglement among atoms and fields in two isolated JCM is done. 

In this paper, we consider the dynamics of the Bloch vector of a two-level atom interacting with a single mode 
boson according to the JCM, where the initial state of the boson is in thermal equilibrium. The time evolution of the 
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Bloch vector S(t) seems to be in a state of complete disorder because of the thermal distribution of the initial state 
of the boson held. The trajectory of S(t) oscillates hard and seems to wander without a purpose. We try to find a 
property that characterizes this confused movement of S(t). 

If S x (0) ^ or S y (0) ^ 0, the time evolution of the Bloch vector S(t) draws a trajectory in the two- or three- 
dimensional space. We observe that this trajectory does not form a closed path. To remove the fast frequency 
oscillation from the trajectory, we take the time-average of the S(t) as (S) = limT->.oo(l/r) f Q T dt S(t). We find that 
(S x ) and (S y ) are equal to zero V/3, w and u>q, where (3 is an inverse of the temperature for the initial thermal state of 
the field, lo is an angular frequency of the field and htoo is an energy gap of the two-level atom. (S z ) takes a certain 
function of /?, u> and ujq. 

We take the histogram of samples {S z (nAt)\n = 0,1, ...,N} for small At and large N. We understand that it 
represents an absolute value of a derivative of the inverse function of S z (t). (When the inverse function of y = S z (t) 
is a multi- valued function, the histogram represents a summation of the absolute values of its derivatives \(d/dy)S~ 1 \ 
at points whose real parts are equal to y on the Riemann surface, where S^ 1 is the inverse function of y = S z (t).) 
We approximate obtained histograms at a high temperature by the probability function of the normal distribution. 
We examine the dependence of the variance of the samples on the temperature (3. 

We examine the time evolution of the entanglement between the atom and the field in our model by estimating the 
entanglement of formation of the density matrix which lies on the 2x2 dimensional projected subspace of the atom 
and the field. Because this projection is an local operation, the entanglement of formation computed in the reduced 
2x2 dimensional subspace gives the lower bound of the entanglement of formation of the whole system (the atom 
and the field). 

This paper is organized as follows: In Sec. [Til we derive an equation which governs the time evolution of the Bloch 
vector S{t). In Sec. IIII1 we examine the trajectory of S(t) and observe that it does not form a closed path. In 
Sec. HVl we take the time-average of S(t) for the case where the field is resonant with the atom. In Sec. [V] we 
take the histogram of data of S z (t) sampled at intervals of At. In Sec. IVI| we take the time-average of S(t) for the 
non-resonant case. In Sec. IVII1 we consider the lower bound of the entanglement of formation between the atom and 
the field. In Sec. IVIII1 we give brief discussions. 



II. THE TIME EVOLUTION OF THE BLOCH VECTOR 



In this section, we give the equation of the time evolution of the Bloch vector of the atom. The Jaynes-Cummings 
model is a system that is described by the following Hamiltonian: 

H = —oj Q a z + hua) a + Tig{a + a + <j-a)), (1) 



where 



a± = -(a x ±ia y ), (2) 



1 _ / —i\ _ ( 1 

1 o J ' ay ~\i or Uz ~ [ o -l 



(3) 



and [a, a^] = 1. The Pauli matrices (cj, i = x,y,z) are operators of the atom and a and a> are operators of the field. 
In this paper, we assume that g is a constant, so that it does not depend on cjq or ui. Let us divide H as follows [HHI]: 

H = h(C 1 +C 2 ), (4) 

d = u(-a z + aU), (5) 

C 2 = g{(J + a + a-a ] ) - -y^ z , (6) 

where Aw = u — ujq . We can confirm 

[C 1 ,C 2 ]=0. (7) 

Because C\ can be diagonalized at ease, we take the following interaction picture. We write a state vector of the 
whole system in the Schrodinger picture as \^>(t)). A state vector in the interaction picture is defined by 

|ft(t)>=exp(iCif)|V(f)>. (8) 
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(We assume |^i(0)) = |^(0)}.) Because of Eq. ([7]), the time evolution of |^>i(t)) is given by 

|^(t)> = Z7(t)]^(0)>, (9) 

where 

U(t) = exp(-iC 2 t). (10) 
We define the density operator of the initial state of the atom and the boson field as 

Paf(0) =pa(0)®/9f, (11) 

p A (0)= Pa,«(0)|*)aa(jI> ( 12 ) 

i,je{o,i} 

exp(— ftfkjja'a) 

PF 



Tr exp(— fthtja'a) 
(1 - e-l 3 ^) exp(-/3?iwa t a), (13) 



where 



|0>A=(J), |1>A = (J) » ( 14 ) 

and ft is an inverse of the temperature. (The indices A and F imply the atom and the field, respectively.) 
The density operator of the atom in the interaction picture evolves according to 

PA(t) = PA,ij(t)\i} AA {j\, (15) 

i,je{o,i} 

PA,ij(t) = PA,kl(®)Akl,ij(t) 

fc,ze{o,i} 

fori,iG{0,l}, (16) 

A kW (t) = A{i\Tr F [U(t)(\k) AA (l\ ® PF )U\t)]\j) A , (17) 

where Trp means a partial trace over the field. Because p At io(t) — p A ,oi(t)* and p A ,u(t) = 1 — p A ,oo(t), we examine 
PA,oo(t) a- n d p A fii(t) only. We can derive Aki t ij(t) as follows: The unitary evolution operator of the whole system 
U(t) given by Eqs. ([6]) and ([10]) is rewritten as 



U(t) = exp[-it 

OO 

= E 



Aw/2 ga 
3a 1, Aw/2 

(-1)"^ 2 " ( (D + g 2 ) n 



n=0 



(2n)! V ° n 



E 



z(-l)»t 2 " +1 / -Alo(D + g 2 ) n /2 gaD r 



^ (2n+l)\ \ ga){D + g 2 ) n AuD n /2 



"00 Uqx 
U W Uu 



(18) 



where 



= (yf + M, (19) 
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. f sm(t^D + g 2 ) 
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u n = cos(tVD) + -Alo — V JL ' . 
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From Eqs. (HZ]), (HH1), (HU) and j20]), we obtain 

(Aw/2) 2 + g 2 (n + 1) cos 2 {\J D(n + l)t) 



A»,oo(*) = (l-e-^)£ 



D(n + 1) 



-n(3huj 



A u ,oo(t) = (l-e-^^ffl ^ e 



sm 2 (yjb(n)t) 



n=l 

oo 



2„ V _____ „~n/3hu 

D(n) 



sin(y / Z)(n+ l)t) 



n=0 \/_D(r. + l) 



i sin(^/_9(ri)t) 



(20) 



x [cos(J D(n)t) - -Aw V ]e-"^, (21) 

2 V^(n) 



4>i,oo(t) = 4io.oo(t) = A»,ox(*) = ^io,ox (*) = A u , i(t) = 0, (22) 

where 

£(n) = (^) 2 + 5 2 n. (23) 

(The trace over the field is taken by the basis vectors of the photon number states.) 
We introduce the Bloch vector S(t) — (S x (t), S y (t), S z (t)) which is given by 

p A (t) = ^(I+S(t)-a), (24) 

where I is the identity operator and er = {(T Xl a yi cr z ). Because pA.{t) — PA(t) and PA(t) > 0, S(t) is a real vector 
and satisfies \S(t)\ 2 < 1. From Eqs. US]), (^2J) and ([2"lj), we obtain 

S(t) = £ Aw (t)S(0) 

-4 2 i(*) I 5(0)+ I ,0 |. (25) 

o o 4 3 i(i) 

where 




L%l(t) = ImLVoi(t)], 

^AuW = ^00,00 (*) — -All, 00 (*)> 

4 4 _W = A 00j oo(t) + ^ii,oo(t)-l- (26) 



This is the equation of the time evolution of S(t). 
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FIG. 1: The trajectory of S(t) (0 < t < 100) whose initial state and temperature are given by S(0) = (1,0,0) and f3 — 0.5, 
respectively. (We assume Aui — 0.) The horizontal and vertical lines represent S x and S x , respectively. They are dimensionless 
quantities. In this case, the trajectory lies on the xz-plane. In the numerical calculation of Li(t) and Li(t) defined in Eqs. (I28p 
and (|30[) . the summations of the index n is carried out up to 500. 



III. THE TRAJECTORY OF THE BLOCH VECTOR 



In this section, we observe the trajectory of the time evolution of S{t). To simplify the discussion, we concentrate 
on the case where the field is resonant with the atom, that is, Aw — 0. Furthermore, we regard ujq as a constant. 
Thus, the model has two variables, t and (3. We replace \g\t with t and (3T%lo(— fthuio) with /3, where we assume g ^ 
and wo ^ 0. They imply that the time t is in units of and the inverse of the temperature (3 is in units of (ftwo) -1 . 
From the above assumptions and Eqs. (|2ip . ([23]) . ([25} and (|26} . we obtain the following equation: 



S(t) = C(t)S(0) 

Li(t) 

L x {t) |S(0)+| I, (27) 







I 3(0) + 









L 3 (t) 
where 

oo 

Li(t) = (l-e-' 3 )^cos(^^n^)cos(^^;^)e~ ,l ' 3 , (28) 



n=0 



1 P 2/3 _ 1 00 

L 3 (t) = ( l- e ^) + _^^cos(2^)e-"' 3 , (29) 



2 V ; 2e^ 

n=l 



L 4 (t) = - - (1 - e"*) + 1 J2 cos(2V^t)e-^. (30) 



(We note Im[j4oi,oi(t)] = for Aw = 0. In the algebraic form of Im[.Aoi.oi(£)], summations of 

cos(V«- + I|ff|£) sm(y/n\g\t)e~ nl3hul / '^Jn and cos(^/n\g\t) sin(V?i + l\g\t)e~ nl3hu} / ^Jn + 1 for the index n from to oo 
appear, where sin(y / n|g|i) / \/n\ n=0 = \g\t. However, they converge to finite values. Similar discussions are explained 
in Sec. II VI ) Looking at Eq. ([27]) . we understand that the trajectory (a curve in the three-dimensional space of S 
parameterized by i) is uniquely determined by the initial state S(0) for given /3. 

Figure Q] shows a trajectory of S(t) whose initial state and temperature are given by 5(0) = (1,0,0) and (3 = 0.5, 
respectively. In this case, S y (t) — for t > and the curve of S(t) lies on the xz-plane. The trajectory of 
Fig. Q] oscillates hard and seems to be in complete disorder. This is because S(t) is a tuple of superpositions of 
cos(^/nT^t) cos(^ht)e- nf} and cos(2^) e - Tl/3 for n = 0, 1, 2, ... as shown in Eqs. ([58} and p)}. 

If the operator C(t) defined in Eq. ([2"7| satisfies the relation £(£2) = £(*2 — t\)£(ti) Vii,i2, the following thing 
can happen. If the trajectory reaches a point (at ia) where it has passed before [at ti(< £2)], it forms a closed 
path and S(t) moves along this closed path for t > t%- (See Fig. O) However, this thing does not happen because 



6 




FIG. 2: The trajectory of S(t) whose initial state is given by S(0) = (1, 0, 0). (We assume Au = 0.) The horizontal and vertical 
lines represent S x and S z , respectively. They are dimensionless quantities. If Citz) = C(t2 — t\)C(t\) Vti,i2, the trajectory 
forms a closed loop. If the trajectory passes a point P at t\ and reaches P at t2(> ii) again, S(t) moves along the closed loop 
for t > t2- This is because the evolution of S(t) is determined by C(t — t2)S(t2). Actually, Cfa) = C(t2 — ii)£(£i) does not 
hold in general and the trajectory does not form the closed path. 
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FIG. 3: The trajectory of S(t) (0 < t < 6) whose initial state and temperature are given by S(0) = (1,0,0) and j3 = 0.5, 
respectively. (We assume Aui — 0.) The horizontal and vertical lines represent S x and S z , respectively. They are dimensionless 
quantities. We can observe that the trajectory intersects itself and does not form a closed loop. In the numerical calculation 
of L\(t) and Li{t) defined in Eqs. (|28[) and (|30|) . the summations of the index n is carried out up to 500. 



C(t 2 ) = C(t2 — t\)£(t\) does not hold in general. In fact, we can observe an example where the trajectory intersects 
itself and does not form a closed path in Fig. [3] 

If S x (0) 7^ or S y (0) 7^ 0, S(t) draws a trajectory in the two- or three-dimensional space. From the above 
considerations, we understand that this trajectory does not form a closed loop in general. 



IV. THE TIME-AVERAGE OF THE BLOCH VECTOR FOR Auj = 



In this section, we take the time-average of the Bloch vector for the case where the field is resonant with the atom 
(Au = 0). As shown in Fig.[lJ the evolution of S(t) is in complete disorder. To remove the fast frequency oscillation, 
we consider the time- average of S(t) as 



1 " T 



T lim - / dt S(t). (31) 
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From Eqs. (27]), (281 (29J) and (301 we obtain 








I S(0) + 






V (La) J 



(32) 



= (1 - e- fJ ) lim — [Vn + 1 sin( Vn + IT) cos(V^T) 

-V^cos(V^TTT)sin(V^T)]e^™' 3 , (33) 

e 2/3 - 1 1 
e _/3 ) H 3— lim — 



n=l 



/r \ X m , (e^- 1 ) 2 ,. l^sin(2V^T) _ n0 



n=l 



Let us take the limit of T — > oo in Eqs. (f33|) . (|34|) and (35]). First, we evaluate (L\), We can obtain the following 
relation about (L%) from Eq. 



-Vx < (Lx) < V u (36) 

where 



Vx = (1-e-l 3 ) lim Iy(V^TT+V^)e 



n=0 



T 1 ™ 00 ^E^ e " n/3 (^°)- ( 37 ) 



n=0 



From Eq. (|37p. we can derive the following relation: 



"X^^rE"^ (38) 

Because 



n=0 



we obtain 



00 ^ 

n=0 ^ J 

^< lim-^. (40) 

Thus, if /3 7^ 0, Vi = and (Li) = 0. Moreover, we can derive (Li) = for /3 = from Eq. (|3"3"]) in direct. Hence, we 
obtain (£1) = V/3. 

Second, we evaluate (L 3 ). We can obtain the following relation about (L 3 ) from Eq. (MJ): 

-,<5!zl,4f«-<,, (41) 

n— 1 v 

where 

^ = ~TJ" 1™ ™ E °)' ( 42 ) 

n— 1 v 

Moreover, we can derive the following relation: 

*<T^J!5. < 43) 

n— 1 
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Because XmLi e — V( e/3 — 1); we obtain 



Thus, we obtain 



By a similar derivation, we can obtain 



1 e 13 + 1 
V 3 <lim ^-j^r = o. 

T—>oc 1 Ae p 



(L 3 ) = -(l-e-P). 



(44) 



(45) 



(46) 



V. THE HISTOGRAM OF S z (t) SAMPLED AT INTERVALS OF At 

In this section, we think about the histogram of {S z (nAt)\n = 0, 1, N} for small At and large N. To simplify 
the discussion, we assume Aw = and 5(0) = (0, 0, 0). The time evolution of Sit) is described as 5(i) = (0, 0, L±{t)), 
where L±[t) is given by Eq. l[5D]l. 

We investigate the time evolution of L±it) by the following way. Fixing (3 at a certain value and defining a small 
time interval At, we collect (JV + 1) samples of L^nAt) for n = 0,1,..., AT, where iV is large enough. Next, we 
make a histogram of these (N + 1) samples {£4(0), L^At), L^NAt)}. We adjust the class interval of bins of the 
histogram, so that the line graph of the histogram approaches a smooth curve. 

This histogram represents an absolute value of a derivative of the inverse function of L^t). As shown in Fig.[H the 
probability that there is a sample L^(n At) in the range of the bin of the histogram [y, y + Ay) is proportional to 



Ll 1 {y + Ay)~L- A 1 {y) 



Ay 



(47) 



where t = L i 1 (y) is an inverse function of y = L^{t) and Ay is a class interval of the bin. (Here, we assume L 4 1 {y) 
is not a multi- valued function.) If we take the limit of small Ay, this probability reaches 



oc 



lim 

Ay-,0 



Ty^ 



L^(y + Ay)-Lf(y) 



Ay 



(48) 



If L7 1 (y) is a multi-valued function, the probability is proportional to a summation of \(d/dy)L^ 1 | at points whose 
real parts are equal to y on the Riemann surface. 

Let us evaluate the histogram of L±{t) in the low temperature limit. In the low temperature limit, we can describe 
L 4 (i) as 



lim L 4 (t) = — (1 - cos2£). 

/3^oo 2 



(49) 



Thus, the histogram is proportional to 



dt 
dy 



v/i-(2y + i) 2 

for -1 < y < 0. 



(50) 



Next, we consider the histogram of L^lt) for the high temperature. If /3 <§; 1, we can expect that L^(t) varies at 
random around a mean with a certain variance. Hence, we can approximate the histogram by the probability function 
of the normal distribution, 



1 



27TCT 



(51) 
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y+Ay 






\. y=L 4 (t) 



FIG. 4: The graph of y = L^it), The horizontal and vertical lines represent t and y (or La), respectively, t is in units of \g\ 
and y (or La) is a dimensionless quantity, t = L^ 1 (y) is an inverse function of y = L^it). 



where fi and a are a mean and a variance of samples, that is, 

AT 

1 

A 4 = 



2V 



1 w 

— ^L 4 (nAt), (52) 



n=0 



and 



1 N 

n=0 

If AT At is large enough, /i is nearly equal to (£4). [From Eq. (f50|) , we can show —1 < L 4 (t) < 0. Thus, we can expect 
that Li{t) varies at random in the range of [—1, 0]. However, Eq. (|5Tj) is defined on —00 < y < 00. Here, we neglect 
this inconsistency.] 

In Figs. 03 [Hand El we show the histogram of L±(t) for /3 = 10, 1 and 0.01, respectively. In Fig. [51 Eq. (|50p fits 
for the histogram well. In Fig. [71 Eq. ([5T|) fits for the histogram well, too. [In Fig. [7J the shape of the histogram is 
not symmetrical. A slope of the left side is steeper than a slope of the right side. The author cannot find a physical 
meaning of this observation.] The histogram of Fig. [S] seems to be an intermediate shape of Figs. [5] and [7J 

In Fig. [51 we plot the variance a 2 of samples against (3 for 0.01 < j3 < 10. For small f3, we can approximate plotted 
points by cr 2 (/3) = ci/3 C2 for 0.01 < (3 < 0.1, where c\ = 0.0516 and c 2 = 2.95. The author cannot find a reason why 
the function of a 2 (/3) for small (3 has such a simple form. 



VI. THE TIME-AVERAGE OF THE BLOCH VECTOR FOR Au / 

In this section, we evaluate the time-average of the Bloch vector for the non-resonant case (Aw ^ 0). Here, we put 
h = 1 (P is in units of fi). We fix g and loq and regard them as constants. Thus, the model has three variables, to, (3 
and t. From Eqs. (|2"Tj). ([2"S"]) . ([2"6"]l and (J3J), we obtain the following relation: 

(L ( H) (L%) 

(S)=| -(4 2 )) |S(0)+| % ), (54) 

(4t 



where 



= (1-e-^) lim I^[ n + ± S in(VJ^ + l)r)co S (VJ>(n)T) 
T ^°° J „=o JD( n + 1) 
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FIG. 5: The histogram of {L4,(nAt)\n = 0, 1, 9999} where = 10, At = 0.05 and the class interval of each bin is equal to 
0.005. The horizontal line represents L4 that is a dimensionless quantity. The vertical line represents the number of samples 
in each bin and it is a dimensionless quantity, as well. A thin line graph represents the histogram of samples {L4(nAt)} . A 
thick curve represents the approximate function o/v/l — (2y + l) 2 where a = 35.5. In the figure, the thick curve is lying on 
the thin line graph and we can hardly distinguish between them. In the numerical calculation of L^{t) defined in Eq. (|30[1 . the 
summation of the index n is carried out up to 1000. 
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FIG. 6: The histogram of {L4,(nAt)\n = 0, 1, ...,9999} where (3 = 1, At = 0.05 and the class interval of each bin is equal to 
0.0035. The horizontal line represents L4 that is a dimensionless quantity. The vertical line represents the number of samples 
in each bin and it is a dimensionless quantity, as well. In the numerical calculation of Li(t) defined in Eq. (|30|l . the summation 
of the index n is carried out up to 1000. 



<4t 



Alo 



D{n) 



-. cos(y/ D(n + 1)T) sin( \J D(n)T)] e - nf3uJ , 



<l-e-^) lim -V 



1 ~ sia(JD(n + 1)T) sm( J D(n)T) 



= (l-e 



+ 



-pus 



n=0 
\'2 



D(n + l)D(n) 



-n(3uj 



D( n y 



DM 

< 2 - 1 ,. 1 g 2 m sin(2 v / J D(m)T) 



UJ 



1 . — i 

. lim 1 y 4 

m—1 



-m(3uj 



=1 D ( m ) y/E(m) 
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FIG. 7: The histogram of \L±{nk£)\n = 0, 1, 9999} where (3 = 0.01, At = 0.05 and the class interval of each bin is equal 
to 1.0 x 10 -5 . The horizontal line represents L4 that is a dimensionless quantity. The vertical line represents the number 
of samples in each bin and it is a dimensionless quantity, as well. A thin line graph represents the histogram of samples 



{L4(nAt)}. A thick curve represents the approximate function [a/(\/2na)] exp 

[~(y ~ M) /(2ff )] where a = 163, \i = -0.00497 



and a 2 = 6.31 x 10~ 8 . (We note (£4} = 
of the index n is carried out up to 1000. 



-0.00498.) In the numerical calculation of L4,(t) defined in Eq. (|30p . the summation 
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FIG. 8: Plots of the variance a 2 of samples {L4(nAt)} against (5 for 0.01 < /3 < 10. Black circles represent plotted data. For 
each black circle, ten thousands samples are taken (n = 0, 1, 9999) and we put At = 0.05. The horizontal line represents /3 
that is in units of (Hljo) -1 ■ The vertical line represents a 2 that is a dimensionless quantity. In both horizontal and vertical 
lines, ticks are put in the logarithmic scale. In the range of 0.01 < [3 < 0.1, plots can be approximated by <r 2 (/3) = cif3 C2 where 
ci = 0.0516 and C2 = 2.95. This approximate function is shown as a line graph in the figure. 



-n(3uj 



D(n) 2^ 'D(n) 

(e^ -l) 2 , 1 ^ g 2 m sin(2 J D(m)T) 
lim — > — 
t — >^ T ^ n, 



4 e /3a 



e —m0u _ ^ 



Using the method shown in Sec. IIV1 we can take the limit T — > 00 in Eq. (|55p and we obtain 



(55) 



l L f u ) = (l-e-^)y[e^^^t + -(e^-l)£^]e- n ^, 
X Aul K D(n) 2 V J D{n) S 



n=l 
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(L<&) = (1 - e-*)Yy»££®- + + l)^r]e-^ - 1. (56) 



Looking at Eq. (|56|) . we obtain 



I™ (ii 3 i) = J™(4 4 i) = 0, (57) 



. ( 3) (Aw) 2 + 2. 9 



and 



The above results imply lim^o^) = VS(0) and 

■ to(S )M°.o, |(A " ) ;r?: i f; ( ° ) " 2a2 ) 

V5(0). (60) 

From Eq. (|60|) . we can say the following: If we start from 5(0) = [the complete mixed state Pa(0) = (1/2) J] in the 
low temperature limit, we can expect to obtain a slightly purified state of S = (0, 0, —2g 2 /[(Auj) 2 + 4g 2 ]) on average 
after the enough time evolution. 



VII. ENTANGLEMENT BETWEEN THE ATOM AND THE THERMAL FIELD 



In this section, we consider the evolution of the lower bound of the entanglement of formation between the atom and 
the thermal field in the JCM. In our model defined in Sec. HH the initial state of the atom and the field is separable. 
However, because of the Jaynes-Cummings interaction, we can expect that the entanglement is generated between 
the atom and the field during the time evolution and their bipartite state becomes inseparable. Such a mechanism 
of entanglement generation is discussed in Refs. fl2l. KL3I1 . as well. Recently, many researchers have been regarded the 
JCM as a source of the entanglement [3, OJ, EMS, El S9, M, H3, H| . 

Let us pursue the time evolution of the entanglement between the atom and the field in our model. To simplify 
the discussion, we assume Auj = and g > 0. Furthermore, we assume the atom initially to be in a pure state of 
(1/\/2)(|0)a + |1)a), which implies the Bloch vector S(0) — (1,0,0) and the density matrix 

Pa(0) = \ ( I I ) . (61) 
The time evolution of the whole state is described as 



Paf(0) = pa(0) ® Pf 
— > pafW - U(t)[ PA (0) ® p F ]U\t), (62) 

where pa(0) is given by Eq. (l6"T|) . p-p is given by Eq. (fl3|) . and U(t) is given by Eq. (fT8|) . 

Here, we arc interested in studying the entanglement for the mixed state of the bipartite system AF (the atom 
and the field). [Because Pf given by Eq. (fT5|) is a mixed state, paf(^) in Eq. (f6"2")l is also a mixed state in general.] 
Moreover, although the dimension of the system A is finite (the two-dimensional system), the dimension of the system 
F is infinite. Entanglement for such a system is difficult to define. 

Some measures of entanglement are proposed at present, for example, the relative entropy of entanglement, en- 
tanglement of formation, and so on. However, because analytical methods are not found for these measures of the 
entanglement in general, it is difficult to compute the value of the entanglement for an arbitrary bipartite state. 
However, exceptionally, an explicit formula of the entanglement of formation for an arbitrary 2x2 state is obtained 
[24l ]. Thus, as the measure of the entanglement, we choose the entanglement of formation. 
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To investigate the entanglement of the bipartite mixed state Paf(£) given by Eq. (|62]) , we take the following method. 
To reduce the dimension of the system F (the field) from an infinite number to a finite number, we project the entire 
state of the atom and the field pAF(t) onto a subspace whose dimension is given by 2 x 2 as 



R AF (t) = (|0) FF (0| 



;ff 



;i|)PAP(i)(io) 



FF \ 



)FF\ 



(63) 



where {|0)f, |1)f} are the photon number states. (We consider the subspace spanned by the basis vectors {|i)a^ |j)f : 
i, j G {0, 1}}.) Because this operation (the projection) is carried out only in the system F locally, the entanglement 
never increases. Hence, the entanglement of -Raf(^) is a lower bound of the entanglement of paf(0. 

After slightly long calculation, we can obtain Raf^) m the form of a matrix which is represented in the basis 
vectors {|0) A |0) F , |0) a |1)f, |1>a|0)f, |1)a|1)f}, 



R AF (t) = -(1 



' r oo,oo r oo,oi r oo,io r oo,n \ 

r oo,oi r oi,oi roi.io r i,ii 

r oo,io r oi,io r io,io rio.n 

\ r 00,ll r 01,ll r 10,ll r n,n / 



(64) 



where 



?"00 


00 


= cos 2 t + e 13 sin 2 t, 


roo 


01 


= ie~P sini cos y/2t, 


roo 


10 


— cost, 


roo 


n 


= — i sin t cos t ~\- ic~ sin t cos t, 


roi 


01 


= e _/3 cos 2 \[2t + e~ 2/3 sin 2 \/2t, 


roi 


10 


= o, 


roi 


11 


= cost cos V%t, 


rio 


10 


= 1, 


rio 


11 


= — isint, 


rn 


11 


= sin 2 t + e~P cos 2 1, 



and we replace phuj with f3 and gt with t. (Here, we assume g > 0. 
Here, for convenience, we rewrite -Raf(£) as follows: 



#AF(*) = PAF{t)(TAF(t), 



where 



and 



Paf(*) = Tr[R AF (t)} 
= (1/4)(1 



)[4 + 3e- /5 + e-^ + e- p (l 



2/3 



J )cos2V2t], 



i^AFft) 

Tr[fl AF (t)]' 



(65) 



(66) 



(67) 



(68) 



[caf (i) is a normalized density matrix of Raf^), that is, Tr[<7AF(t)] = 1.] 

The original definition of the entanglement of formation of an arbitrary density matrix pabi which is a bipartite state 
of systems A and B, is as follows. Let us suppose pab is shared by Alice and Bob. And suppose that, asymptotically 
as n — > oo, Alice and Bob can prepare (pab)®" from k Bell pairs using local operations and classical communication. 
The entanglement of formation of pab is given by 

E(pab) = lim (69) 

n — >oc Ji 

where /c m in is the minimum of k for given n [25| . 

An explicit formula of the entanglement of formation of a 2 x 2 dimensional (normalized) bipartite density matrix 
Pab is given as follows. First, we compute the concurrence of pab from which we can calculate the entanglement of 
formation of pab- We define a matrix pab as 



Pab = (o-y ®<y y )p* AB (a v 



(70) 
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FIG. 9: The variation of E(Rpj?(t)) against t G [0, 2ir] with fixed f3 — 10, 2 and 1. The horizontal and vertical lines represent 
t (time) and E (entanglement of formation), respectively. The time t is in units of g -1 . The entanglement of formation E is a 
dimensionless quantity. A thick curve, a thin curve, and a dashed thin curve represent the variations of E(Rp l p(t)) for f3 — 10, 
2, and 1, respectively. 



where p AB is the complex conjugate of pab that is represented in the basis vectors of {|*)a|j)b : hj € {0, 1}}- We 
write the eigenvalues of PabPab as Ai, A2, A3 and A4, where Ai > A2 > A3 > A-i(> 0). The concurrence C(/?ab) is 
given by 



C(pab) = max{0, y/Xi - \/A 2 - \/A 3 - \/A 4 }. (71) 
The entanglement of formation of pab is given by £ (C(pab)), where 



E(P) = 2 g2 2 2 g2 2 ' ^ ' 

Here, we want to estimate the entanglement of formation of i?AF (t) = Paf(£)caf(£), which is not a normal- 
ized density matrix. The entanglement of formation of i?AF(£), E(RAv{t)), should be written as E(RAF(t)) = 
Paf(£)£(C(oaf (£)))■ Because the analytical form of <7af(*) obtained from Eqs. ([6"4"|) , ([55)) and ([fJ5|) is very com- 
plicated, it is difficult to obtain C(oaf(£)) m an explicit formula. Thus, we estimate C(ctaf(£)) [and £ , (i?AF(^))] 
numerically. The variations of E(B,AF(t)) against t G [0, 27r] with fixed /3 = 10, 2 and 1 are shown in Fig. [5] 

In Fig. [SI 2?(-Raf(*)) against £ seems to vary periodically. Because we concentrate on the photon number states 
|0)f and |1)f only, the effect of the thermal distribution of the field is neglected. Looking at Fig. [HI we notice that if 
P takes a large value (at a low temperature), a certain amount of the entanglement arises between the atom and the 
field. As (3 becomes smaller (as the temperature becomes higher), the amplitude of oscillation of E(B,AF(t)) becomes 
smaller. However, this observation does not imply that the entanglement of formation of the entire system (the atom 
and the field) becomes smaller as the temperature becomes higher, because E(R,AF(t)) is just the lower bound of 
the entanglement of formation of paf(£)- From these results, we may regard the thermal JCM as the source of the 
entanglement. 



VIII. DISCUSSIONS 



In this paper, we examine the dynamics of the Bloch vector of the two-level atom in the thermal Jaynes- 
Cummings model (JCM). In the evolution of the Bloch vector, for example, if Au> = 0, infinite summation of 
cos(v / ^ + 1£) cos(v / n£)e~™' 3 and cos(2y / n£)e - ™' 3 for n = 0,1,2,... appears and this makes it difficult to treat the 
problem exactly. 

In our model, the thermal effects are introduced only in the initial state of the boson field. To discuss the thermo- 
dynamics of the JCM strictly, we have to think the grand partition function of the whole system (the atom and the 
boson field) and pursue its non-equilibrium time evolution. Although the JCM has been studied by many researchers, 
understanding about the thermal JCM seems not to be enough. 
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In this paper, we try to obtain a global property that characterizes the confused behavior of the Bloch vector. We 
observe the trajectory of the Bloch vector and take its time-average. We take the histogram of the z components 
of the Bloch vector sampled at intervals of At. However, the author wonders whether these results are good global 
aspects of the trajectory that is in complete disorder. 

Recently, entanglement generation during the evolution of JCM has been studied from the viewpoint of the quantum 



information thcor; 
also discussed 



eory I!. 

MM 



F2JJ13 
2111 



IH Kdi . [I7L Hil]. On the other hand, the entanglement sudden death in the JCM is 
231 1. These matters are explained in Sec. HI In Sec. IVII1 we consider the evolution of 



the entanglement of formation between the atom and the thermal field. The thermal JCM may become an important 
source of the entanglement in future. 
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